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ABSTRACT 

The PHOTOS Monte Carlo is often used for simulation of QED effects in decay of intermedi- 
ate particles and resonances. Momenta are generated in such a way that samples of events cover 
the whole bremsstrahlung phase space. With the help of selection cuts, experimental acceptance 
i-C^ can be then taken into account. 

The program is based on an exact multiphoton phase space. Crude matrix element is ob- 
tained by iteration of a universal multidimensional kernel. It ensures exact distribution in the 
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^ soft photon region. Algorithm is compatible with exclusive exponentiation. To evaluate the 
^ program's precision, it is necessary to control the kernel with the help of perturbative results. If 
■ availabe, kernel is constructed from the exact first order matrix element. This ensures that all 
\^ ' terms necessary for non-leading logarithms are taken into account. In the present paper we will 
O focus on the W ^ l\ and 7* — )• tt+tt^ decays. The Born level cross sections for both processes 
^ approach zero in some points of the phase space. 

A process dependent compensating weight is constructed to incorporate the exact matrix 
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element, but is recommended for use in tests only. In the hard photon region, where scalar QED 
is not expected to be reliable, the compensating weight for 7* decay can be large. With respect 
5^ to the total rate, the effect remains at the permille level. It is nonetheless of interest. The terms 
leading to the effect are analogous to some terms appearing in QCD. 

The present paper can be understood either as a contribution to discussion on how to match 
two coUinear emission chains resulting from charged sources in a way compatible with the exact 
and complete phase space, exclusive exponentiation and the first order matrix element of QED 
(scalar QED), or as the practical study of predictions for accelerator experiments. 
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1 Introduction 



One of the crucial goals of any high energy physics experiments is the comparison between 
results of new measurements and predictions obtained from theory. If agreement is obtained, 
then the validity domain for the theory is extended. Discrepancy can be attributed to the so 
called new physics. This scheme is in principle rather simple, but in practice, it is involved. For 
LEP experiments, enormous effort for such a program was documented in |[I1I3- It was nec- 
essary because for very precise scattering experiments one needs to study radiative corrections 
simultaneously with detector acceptance. As a consequence, it was possible to confirm experi- 
mentally that the Standard Model was indeed a field theory of elementary particle interactions. 
Quantum effects could not be omitted, they had to be included in calculations. The importance 
of this achievement was confirmed by the 1999 year Nobel Prize attributed to 't Hooft and 
Veltman. Also the 2008 Nobel Prize for the mechanism of quark flavour mixing [3J required 
precise measurements and comparison with data. These experiments. Belle and BaBar, located 
respectively in Japan and USA required good control of radiative corrections too [|4l|51. 

Because of the nature of accelerator experiments, it is generally believed that the Monte 
Carlo technique is the only suitable tool for the precision comparison of theory with experi- 
mental data |T|: effects of detector acceptance can be merged into theoretical predictions by 
simple rejection of some of the generated events. Theoretical effects of different nature can be 
taken into account, in particular radiative corrections. A multitude of Monte Carlo programs 
were developed in context of QED [|6l|7l and QCD [Hill. The references can serve as examples. 

Such Monte Carlo programs must rely on results obtained from perturbative methods. In 
the case of QED, exponentiation [lOJ is useful. Exponentiation is a long established and rigor- 
ous scheme of reorganization of perturbative expansions. It was found fTT| that Monte Carlo 
programs can be developed using it as a basis. Significant theoretical effort was nonetheless 
necessary. It required not only explicit calculation of exact fixed order cross sections, but also 
to separate them into appropriate parts, at the cross section or the spin amplitude level, to finally 
match results of fixed order calculations with coherent exclusive exponentiation (CEEX) [|T2l| 
and implement it into computer programs. 

QED predicts distributions which are strongly peaked in phase space. They may vary by 
more than 10 orders of magnitude. This and the complex structure of infrared singularity can- 
cellations, poses a challenge in Monte Carlo method. Appropriate choice of the crude distribu- 
tion over the phase space must be found. In the case of multi-photon radiation, a particularly 
elegant method was found [fT2| . Thanks to conformal symmetry, it was possible to construct 
a crude distribution which was actually exact from the phase space point of view. All simpli- 
fications were localized in an approximated matrix element. An approximated matrix element 
consisting of the Bom amplitude multiplied by the so-called soft factors was used at a first step. 
Any further improvements could then be easily achieved with a correcting weight. The weight 
is the ratio of distributions calculated from an available matrix element obtained perturbatively 
to a given fixed order and the one used in first step of the generation. 

The case of QCD is by far more complex, but the general principle is similar. One constructs 
a simplified matrix element (and approximated phase space) as a basis for the parton shower 
algorithms. Such a solution is limited to leading logarithms [fT3] - [T7l . Phase space organizations 
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based on so-called orderings are often used HHIHl . Improvements beyond leading approxima- 
tions are possible and widely used II20[|21L but have technical difficulties, for example through 
the appearance of negative weight events. 

Another difficulty of QCD is the necessity to use parameterization of a., at low where the 
standard parameterization a., = (2^) becomes unreliableQ. At small scales non-perturbative 
aspects of QCD can dominate. Phenomena like underlying event [|24l or hadronization [|25ll 
lead to further complications. This is of course on top of hadron structure functions obtained 
from experiments f2F|. 

We will not elaborate in our paper on these topics, we think however, that the methods ap- 
plied in this paper may provide a useful hint. First results [27] are encouraging. Our paper is 
devoted to reliability of the PHOTOS Monte Carlo. PHOTOS is a Monte Carlo I128IJ29I for the QED 
bremsstrahlung in decays. Its structure is similar to the algorithms for QED exclusive exponen- 
tiation; the parameterization of its phase space is exact, but its algorithm is iterative!!. Conformal 
symmetry is not used. This is advantageous, terms responsable for leading logarithms of decay 
product are reproduced to all orders. 

As in the cases mentioned before effort in understanding results of exact perturbative cal- 
culations was necessary for the construction of PHOTOS. Cross section level distributions, Refs. 
|[30l[3T1l . were used. Later, thanks to experience gained in KKMC project [[T21I321 . spin ampli- 
tudes were found to be helpful. In particular, results of Refs. |[33l - [35l were used. They were 
essential for design and tests of the program, in particular for the choice of single emission ker- 
nels. Thanks to these works, interference of consecutive emissions from a charged line as well 
as interference of emissions from distinct charged lines was properly taken into account, with- 
out any need to divide the phase space into differently treated sectors. That is also why there 
was no need to separate the photon emission phase space into regions where either a shower or 
a fixed order hard matrix element is used. 

Refs. [l36] - [38l were devoted to numerical tests, but also a better explanation of theoretical 
foundation of the program was given there. It may be worth mentioning that for many years 
the program's precision was of no interest and such explanations were delayed to the present 
decade. 

The best detailed description of the phase space parameterization as used in PHOTOS and the 
explanation that it is actually exact, is given in Refs. Il37l[38l . However, it is not different from 
what was already explained in 1128 [|29l . 

The precision of the program is significantly improved with respect to its early versions. As 
it was shown in [|36[|3711 . even if the incomplete first order matrix element is used in PHOTOS, its 
results agree much better with KKMC using second order matrix element exclusive exponentia- 
tion [|32ll than with KKMC, using matrix element restricted to first order and exponentiation. To 
quantify this statement, the method described in B9IJ40II was used. PHOTOS was found to exploit 
result of perturbative calculation quite well, but it cannot be a substitute of such calculations. 

The main goal of the present paper is to study spin amplitudes for construction of the process 

' Discussion of this problem can be found in 022II23I . 

^That is why it is similar to solutions used in QCD parton showers, but no phase space ordering of any sort is 
applied and of course most of the difficulties present in QCD are absent as well. One should also stress differences; 
iterated, single emission kernel simultaneously feature all emission sources. 
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dependent weight, as in B7II38II . but for 7* — tt+tt^ and — /^v. For that purpose, the matrix 
element needs to be studied in great detail. Its gauge invariant part^ need to be identified and 
with their help relations with amplitudes of lower orders have to be found. This second aspect 
is important and is closely related to properties used in defining factorization schemes, see 

eg. mm. 

Let us point out that in this paper we will not discuss spin amplitudes from the perspective 
of matching consecutive emissions from the same charged line. Such studies were performed 
earlier ||29ll and for other decays; these studies required double emission QED amplitudes 041 
[35l . We will focus on single photon emission and matching the emissions from two charged 
lines in 7* — )• 71+71". The analysis of the spin amplitudes and tests for the algorithm in the case 
of Z decay into pair of charged fermion was given earlier, in ref. Il37l . The scalar particle decay 
into a pair of fermions was covered in Il43l and the decay of a spinless particle into a pair of 
scalars was studied in [38J. It seems that the algorithm works better (correction weights are 
less important) when initial state is spinles^. The case of W decay was covered in [,441 . though 
some approximations were used and the decay requires to be revisited. 

The two processes are not only of the technical interest, they provide examples for studies 
of Lorentz and gauge group properties of spin amplitudes and cross sections. The y* — )■ 7r+7r^ 
decay is well measured. It is important to improve theoretical uncertainty of PHOTOS for this 
decay, because of its relevance to establishing <Xqed{Mz) and to phenomenologvof g — 2. From 
that perspective, the validity of our study is limited by validity of scalar QEEO. The W ^ l\ 
decay is of interest for precision measurement of W mass and width, at LHC for example. 

Our paper is organized as follows. In Section 2 we present the scalar QED spin amplitudes 
for the process e+e^ — )• 7* — )■ 71+71^(7). It will be shown that the spin amplitudes can be sepa- 
rated into two gauge invariant parts. Section 3 covers further discussion of the amplitudes, and 
the formulae for the cross section used to obtain numerical results. A separation into eikonal 
part and remaining parts is also presented. Section 4 is devoted to the numerical results ob- 
tained with the help of MC-TESTER [|39l|40l. Different options of separating non-leading effects 
are demonstrated. Section 5 is devoted to the discussion of further tests, where distributions 
sensitive to the beam direction will be used as well. Similarities and differences with respect 
to the previous case will be underlined. Section 6 summarize the paper. Spin amplitudes for 
— )• /+v are given in the Appendix. 

2 Amplitudes 

One of the necessary steps in the development of any Monte Carlo program is analyzing spin 
amplitudes calculated from the theory or phenomenological model under consideration. Fixed 
order analytical results are often not sufficient. Even well known amplitudes have to be revisited 
again to study their structure. It appears to be fruitful to study decompositions of the amplitude 

^ In our case, gauge invariance reduces to independence of longitudinal component of photon polarization. 
^This is particularly interesting from the point of view of future attempts to extent into QCD. 
^ This last constraint is of course common with the projects such as PHOKARA 1451 . PHOTOS will not be better 
or worse from that point of view. 
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into sums of gauge invariant parts, which can be further factorized into gauge invariant terms. In 
particular, to find correspondence with factorization properties of the underlying field theorjl^. 

The spin amplitudes for W — t- Ivy are collected in the Appendix. In principle they are 
straightforward and available already in appropriate form in [44] but let us recall them again to 
clarify possible ambiguities on how emission from W is separated into final state radiation and 
initial state radiation. Let us point also to [|47ll . where spin amplitudes for radiative corrections 
and W decay Monte Carlo are discussed for the first time. This is an important starting point 
for the discussion of radiative corrections necessary for precise measurement of W lineshape. 

In the following, let us concentrate on the process 'fip) ^^(^i)^ (^2)Y(^,e)- Since the 
precision required by experiments is lower in this case than at LEP, we will limit ourselves to 
the discussion of amplitudes for single photon emission. We will not perform detailed analysis 
of virtual corrections. At required level of precision, it is enough to anticipate their size thanks 
to the Kinoshita-Lee-Nauenberg theorem Il48ll49ll . Anyway, scalar QED predictions for our 
process are only partly reliable. 

The spin structure of our process is new with respect to the processes we have already 
studied: Z — )• Z+/^, h — )• Z+/^ or — )• tc+tc^. The spin of the initial state can not be transmitted 
into helicities of the outgoing particles. That is why we can expect different properties of the 
amplitudes. 

If one considers the process e^e^ Tip) ~^ ^^(^i)^ (^2)y(^, e), amplitudes equivalent 
to those given in [|50ll are obtained. The amplitude can be written as M = V^/Zu where = 
v{pi, )y^u{p2, '^2). The pi , Xi , p2, ^2 are momenta and helicities of the incoming electron and 
positron. The define the spin state of the intermediate 7* . 

Let us turn to the virtual photon decay now. Following conventions of [50], the final inter- 
action part of the Born matrix element for such process is 

H^ip.qi.qi) = ^^^(^1 -^2)^. (1) 
Here p = qi+q2- If a photon is present, this part of the amplitude reads: 



P i qi-k q2-k J 



(2) 



which can be re- written to the following form: 



q\-k q2-k 
e^F2n{p^) [k^qi ■ £* - e*^^i ■ k k^q2 ■ £* - £*^q2 ■ k 



+ ^ -^^\. (3) 



P \ qi-k q2-k 

Formally, //o(p,^i,^2) is as in the Bom case, but with p = qi+q2 + k instead of p = qi+q2 
for the virtual photon propagator. Let us note, that the first Born-like term and two other terms 



*The particulary rich case of e+e VeVeY and Monte Carlo implementation of spin amplitudes separated 
into parts is discussed in | 
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in the second line of ^ are separately gauge invariant. The normalization of space-like part of 
Hq is not as at the Born level: |^i — ^2! < The time-like part of Hq drops out when 

the product with Vu is taken. 

Let us now return to Vy. Following conventions of fSTI it reads: 

V^ = 2{\X+\e^^ + iX+e'^ - e^) . (4) 

where X± = Xi ±X2. The vectors satisfy {eaY = 5^. We choose ei to lie in the reaction plane, 
while ^2 = Pi X (^1 — <l2)/\pi X (^1 — 12) I is chosen to be perpendicular to that plane. The ei 
is along incoming electron beam and eo is proportional to pi + p2- The basis vectors can be 
written as 

eo = ( 1 ) ,ei = ( 1 ) ,^2 = ( 1 ) ,63 = ( 1 ) . (5) 

We can drop the term proportional to the electron mass. 

At Born level the second term in the expression dl]) will not contribute because 62 ■ {qi — 
^2) = 0- The complete amplitude is thus: 

2 1 

MBorn = e F2Ti{S)-j=\X^\ei ■ {qi-q2), (6) 

where a/S is the energy of cm.. One can see that the amplitude is proportional to sinGg as it 
should be. Here 0g = Zpi^i is a scattering angle. Squared and summed over initial spin states 
the amplitude yelds: 

V- , ,9. X 8(47ra)2F2 (5) ... 
Y,\MBorn\Hs,T,U) = ^ ^2 {TU-mls) . (7) 



The Mandelstam variables are defined as follows: 



S = 2pi-p2, T = 2pi-qu U = 2pi-q2. (8) 

The amplitude ([3]) for single photon emission can be decomposed into a sum of two gauge 
invariant parts: 

Hf' = Hf + Hfj (9) 



or 



where 



H''^H^,+H^j, (10) 
e^F2n{p^) ^ ,^f qi-£* qi-t y .... 



K = f ^ + ^ V , (12) 

p^ \ \qi-k q2-k J J 



and alternatively 



p^ V q2-k + qi-k/ \ qi-k q2-k ' 



2e^F2n{p^ f ^ 

p2 \q2-k + qi-k' 



Hi, = ' ( - >- - ,. (gi-e^ + g2-e1-e^n ■ d^) 
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3 Cross section and amplitude separation 



Before going into numerical results, let us elaborate on the formulas presented above in more 
details. One can see rather easily that formulas (fTTI) and (fT3l) have a form typical for amplitudes 
of QED exclusive exponentiation Il52ll . that is, Bom factors multiplied by an eikonal factor 

{^yk ~ ^qrk) ' ^^'^ expressions differ in the way how the Bom factor approach the 

genuine Born expression. In both cases the expressions approach the Bom in soft photon limit. 
In case of (fT3l) this property holds for the photon coUinear to 71+ or . This was achieved 
by adding to (fT3l) the term proportional to ^'" ^^'^^^''^ and subtracting it from (fT4l) . As a con- 
sequence the expression in the first bracket of (fTfl . in coUinear configurations will be close 
to ^1 — ^2 =F ^ respectively if qi ■ k <^ q2 ■ k and q2 ■ k <^ q\ ■ k. Thus, it is consistent with LL 
level factorization into Born amplitude and eikonal factor. Generally expressions (fTTI) and (fT3] ) 
differ from a product of Born times eikonal factor only by normalization. This defect is easy 
to correct, and we will return to this point later in this section when discussion of cross section 
will be given. 

Experience with the Z — )■ /+Z^ decay has shown that it is useful not only to rely on spin 
amplitudes, but to collect expressions for amplitudes squared and (partly) averaged over the spin 
degrees of freedom, since it can be useful for future work on matching kernels of consecutive 
emissions. 

If one takes separation ^ for the calculation of two parts of spin amplitudes, after spin 
average, an expression for the cross section based on ^ takes the form: 

£|M|2 = £|M/|2 + J^|M//|2 + 2£M/M;,, (15) 

A,,e A,,e A.,e A,,e 



where 



l^\M,\ = -2e + 



^ 52 \{qi-k)^ {qi-kY {qi-k){q2-k) 

{ {TU' - mis) + {t'U - mis) - {TT' + UU' - SS') ] 
2 .^2 5/ 



8^'^ { tAt, + tAt, - . Z ] (^1 • k) (^2 ■ kl (16) 



52 \{qi-k)2 {q2-kf {qi-k){q2-k) 
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VIM |2 - 0^6 -PfJ'^) f f ml S' 



52 \\iqi-k)^ {qi-kf iqi-k){q2-k) 
{{TU' - mis) + {T'U- mis) + {TT' + UU'- SS') ) - 85 



4 

4 
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— (jU' - mis + ^{TT' + UU' - SS')^ 
(j'U - mis + ^{TT' + UU' - SS')^ | 



, .F^T,{S) ( f ml ml S' 



52 {\{qi-k)2 {q2-kY {qi-k){q2-k) 
(A+5-C-D) + -1^(54(C + D))+-1^(a-1(C + D) 



{^_{TU' -mls) + {T'U -mis)) 
2 



2 



'— (tu' -mlS-^ {TT' + UU' - SS')^ 
— (^-{T'U -mls) + ^{TT' + UU' - SS')^^ 



{q\-kY {qi-kYJ q\-k 



52 V(^l-^)' {qi-kf {qi-k){q2-k) 
+Se''^^{S' + 2ml + S), (17) 



^2 



A + -(C + D) +8e^^|^(5' + 2m^-5). (18) 



q2-k V 2' 7 J 52 

The definitions of terms A,B,C,D,E will be given later in the section. 

If instead of expression Q we use (flOl) the following relation is obtained: 



£|M|2 = £|M,,|2 + £|M„,|2 + 2£m,,m;,,, (19) 



1 



where 



i(^i-fc + ^2-^)2 ^ ^ {qi-k + qi-kf^ ' 



{qi-k + q2-k) 



.2 ^2 



= -{4%a){-. — + - 



{qi-kY {q2-kY {qx-k){q2-k) 



{qi-k + q2-k)^ {qi-k + q2-kY {qi-k + q2-k)^ 

...6FliS)f ml ml S' \ {qi-kf{q2-kf , 

52 \{qx-kf^ {q2-kY {qx-k){q2-k)) {q,-k + q2-kf ^ 



-^e^FliS) 



52 



{q\-k + q2-kY 



{TU' - mis + T'U - mis + TT' + UU' - SS') 



-2S] 



{q\-k + q2-kY 



52 {qi-k + q2-k)^ 



(21) 



Note that Y,x,e Wii'l^ is free of infrared and coUinear divergences. The interference contribution 
is given by the following expression: 



-^e'FliS) 



52 

(Iml + S 



{qi-k + q2-k)^ 
A G2 ■ k 



i2ml + S')^^-S]{TU'-mls) + 



q\ -k 

S / qi-k ^ q2-k 
-Ana 



S^ {T'U -mis) + 
-2ml-S']{TT' + UU'- SS') 



{2ml + S') 



j^(l2-k 



{q,.k + q2-kY\y--'^'qr.k ")^+((^'"^ 



+ S')^-S\B 

qi-k 



S^{qi-k + q2-kY 



{{qi-kf + {q2 -kf)- {2ml + S') {qi'k) {q2 ■ k) 



(22) 
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The Mandelstam variables are defined as follows 



Finally 



S = 2p\-p2, S' = 2q\-q2, 
T = 2pi-qi, T' = 2p2-q2, 
U = 2pi-q2, U' = 2p2-qi. 

A = £|Mfi„,„|2(5,r',[/), 
A, 

B = Y,\MBorn?{SJ,U'), 

X 

c = £|MB„,„|2(5,r,t/), 

A 

D = Y,\MBorn\^{Sj',U'), 
A 

£ = 32(4m)3m2^. 



(23) 
(24) 
(25) 

(26) 
(27) 
(28) 
(29) 

(30) 



Let us point that the complete expression for the amplitude squared is, in comparison to its 
parts, short: 

We should stress that our two separation options (eqs.(fT5l) and (fT9l) ) can have their first terms 
even closer to Bom-times-eikonal-factor form. For that purpose it is enough to adjust normal- 
ization of (fTTI) (or (fT3l) ) to Bom amplitude times eikonal factor. Compensating adjustment to 
(fT2l) (or (fT4l) ) is then necessary. The changes can be performed by numerical manipulation of 
the three contributions to ([T?]) and (fT9l ). The resulting new separation into parts will be dis- 
tinguished by additional prime over its parts. For example Y^x ^^ will be used instead of 

Such a modification is of interest, because if 'Yx^^ l^/P or E?i,e l^/'P is xx^qA alone, then it 
is the expression used in simulation with PHOTOS Monte Carlo and refinement of [38]. In the 
next section, we will perform our numerical investigations with respect to results obtained from 
formulas of ref. Il38ll (which for our present process is just an approximation). 



4 General numerical results 

We have performed our numerical studies for the decaying photon virtualities of 2, 20, 200 and 
2000 GeV. However in this paper we will show only the case of 2 GeV. The other ones confirm 
only that the coUinear logarithms are properly reproduced by the simulation with standard set- 
up of the PHOTOS kernel and would not add anything relevant to our discussion. 
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Let us start with a presentation of the case when the weight for the matrix element is that 
for Y^x ^ l^/P or Y.'x^e 1^/' one can see from fig {!] agreement with matrix element of Il38l is 
excellent over all the phase space. Unfortunately, as it will be discussed in the next section, this 
is true only for the case when distributions are averaged over the orientation of the whole event 
with respect to incoming beams (or spin state of the virtual photon). At this moment let us note 
that as a consequence of strongly varying Born cross section (approaching zero in forward and 
backward direction) the resulting weight distribution from Y,x,e I^/P or E?i,e has a tail. 
We have used special techniques to appropriately adapt Monte Carlo simulation to that. 

If Ll.e \Mi\^ or Y,x,E \Mji\^ is used directly instead of g or Y^x.e normalization 
of Born-like factor is not performed, differences with respect to formulas of [|38l are much 
larger, see respectively figs [21 and [3l In the last case discrepancies are smaller, because the 
normalization is correct in coUinear limit. Finally let us compare results of complete scalar 
QED matrix element with that of [38J, see fig. IH In the high photon energy region there is a 
clear surplus of events with respect to the formulas of [,38,1 . That contribution should not be 
understood as bremsstrahlung, but rather as a genuine process. Anyway in that region of phase 
space scalar QED is not expected to work well. Even though expression (|3TI) looks elegant and 
is short, it needs to be separated into (at the cross section level) longer expressions, where Bom 
times eikonal factor part is explicitly separated. Note the difference between results shown on 
figs [H and |4] is only 0.2 % of the total rate. That is why our detailed discussion is not important 
for numerical conclusions, but for the understanding of the underlying structure of distributions. 
Once the status of approximations used in PHOTOS at single photon radiation is understood, we 
can, as in other processes iterate and to simulate effects of multibremsstrahlung simultaneously 
with the detector effects. As an example we show in fig. [5]results of the single photon emission 
mode, and compare with the one of multiple emission mode. Differences are rather small. This 
may not be the case if selection cuts are present. 

Now, let us consider the decay W — )■ /v(y) . In [|44ll a simple correcting weight was introduced 
into PHOTOS for discrepancies with respect to exact predictions of SANC ll53l . A weight based 
on the exact matrix element is presently available, see Appendix. One can check again how 
good approximation of ref. [44] is. As one can see from fig. [6] the correction weight reproduces 
the result of exact matrix element well. In fig. U\ we show that once the exact matrix element 
is implemented into PHOTOS the agreement with the benchmark calculation is better than the 
statistical error of 10^ events. In contrary to the previously studied y* TT+TT^y case, there 
were no problems of weight distribution tail. 

5 Numerical results using beam direction 

In the previous section we have discussed distributions regarding four-momenta of decay prod- 
ucts only. Agreement between results of PHOTOS using universal kernel and simulations based 
on matrix element was excellent both in case ofW—> ^v(y) and 7* 71+71^(7) decays, even 
though the decaying particle spin effects were not taken into accounl|Z|in the PHOTOS kernel. 

^This is definitely a complication requiring some attention. It is an interesting aspect of the validation of 
PHOTOS, absent in the scalar state Il38l . but present in case of Z decay [37], and it is strongly related to limits of 
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Figure 1: Distributions of invariant masses normalized to center of mass energy and squared (M'^ jS) 
for e^e~ — JI"'"JI~(y) at 2 GeV center of mass energy. Results from PHOTOS with matrix element taken 
from l\38V are given in red (or darker grey) colour If matrix element Y,x.e or Y.Xe I^/'P used (the 
two options are effectively identical) results are given in green colour Logarithmic scale is used, but 
for the ratio (black line) linear scale is used instead. Fraction of events with photons above 50 MeV is 
respectively 4.2279 ± 0.0021 % and 4.2269 ± 0.0021% for the two programs. 




(a) This distribution is identical to the distribution of (b) This distribution coincides with distributions for 
photon energy in the reaction frame as well. squared invariant mass ofn^y pair. 



Figure 2: Distributions of invariant masses normalized to center of mass energy and squared (M^/S) 
for e^e^ — >• 7I^7I^(y) at 2 GeV center of mass energy. Results from PHOTOS with matrix element taken 
from l\38\l are given in red colour If matrix element \Mj\^ is used results are given in green colour 
Logarithmic scale is used, but for the ratio (black line) linear scale is used instead. Fraction of events 
with photons above 50 MeV is respectively 4.2279 ± 0.0021 % and 3.4435 ± 0.0019 % for the two 
programs. 




(a) This distribution is identical to the distribution of (b) It coincides with distributions for squared invariant 
photon energy in the reaction frame as well. mass ofn^y pair. 
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Figure 3: Distributions of invariant masses normalized to center of mass energy and squared (M'^ jS) 
for e^e~ — >■ JI"'"7I~(y) at 2 GeV center of mass energy. Results from PHOTOS with matrix element taken 
from USS^l are given in red colour If matrix element Y,x.e |^/' P used results are given in green colour 
Logarithmic scale is used, but for the ratio (black line) linear scale is used instead. Fraction of events 
with photons above 50 MeV is respectively 4.2279 ± 0.0021 % and 3.8329 ± 0.0020 % for the two 
programs. 




(a) This distribution is identical to the distribution of (b) It coincides with distributions for squared invariant 
photon energy in the reaction frame as well. mass ofn^y pair. 



Figure 4: Distributions of invariant masses normalized to center of mass energy and squared (M'^ /S) 
for e^e^ — >• 7I^7I^(y) at 2 GeV center of mass energy. Results from PHOTOS with matrix element taken 
from 081/ are given in red colour If complete matrix element is used results are given in green colour 
Logarithmic scale is used, but for the ratio (black line) linear scale is used instead. Fraction of events 
with photons above 50 MeV is respectively 4.2279 ± 0.0021 % and 4.4320 ± 0.0021% for the two 
programs. 




(a) This distribution is identical to the distribution of (b) It coincides with distributions for squared invariant 
photon energy in the reaction frame as well. mass ofn^y pair. 
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Figure 5: Distributions of invariant masses normalized to center of mass energy and squared (M"^ jS) 
for e^e~ — >■ TC"''7I^(y) at 2 GeV center of mass energy. Results from PHOTOS with matrix element taken 
from ussy are given in red colour If exponentiation is activated as well results are given in green colour 
Logarithmic scale is used, but for the ratio (black line) linear scale is used instead. Fraction of events 
with at least one photon above 50 MeV is respectively 4.2279 ± 0.0021 % and 4.1377 ± 0.0020% for 
the two cases. 




(a) This distribution is identical to the distribution of (b) It coincides with distributions for squared invariant 
photon energy in the reaction frame as well. mass ofn^y pair. 



As can be seen from plots ISUTTl distributions in variables sensitive to the orientation of 
the W boson spin are affected. The plots [8] and [9] show distributions in the photon momentum 
angle with respect to a spin quantization axis as predicted by SANC and by PHOTOS with the 
standard kernel in transversally and longitudinally polarized W boson decays. The plots \T0\ 
and \TT\ correspond to the muon momentum orientation. The regions of phase space, where 
distributions are sparcely populated and where in fact at Bom level probability density approach 
zero, are becoming moderately overpopulated by PHOTOS (increase of up to 14 % of density 
was found for transversely polarized W boson decay). In most cases, this is probably of no 
practical consequences, nonetheless it requires quantification. Once the exact matrix element is 
implemented into PHOTOS, agreement with the SANC predictions is better than statistical error of 
10^ events, see fig.[T2l 

Similar effects take place for y* — )• TT+TT^y. Even though from fig. [Done could conclude that 
the universal kernel of [38], for arbitrary large samples, is equivalent to the matrix element as 
given by Y^x.e l^/P or J^^.e 1^/' differences appear in distributions sensitive to initial state spin 
orientation, see figs \T3\ and [141 On these plots angular distributions of the photon momentum 
with respect to the beam line are shown. Again, regions of phase space giving zero contribution 
at the Born level are becoming overpopulated if an approximation for the photon radiation 
matrix element is used. From that perspective and for practical reasons one can conclude that the 
L?i,e 1^/' P choice is better than Y^x.e It yelds distributions closer to the ones obtained from 
universal kernel. Then, the remaining part of ([T9l ) represents better correction to implement bulk 

factorization. 
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Figure 6: Ratios of results from PHOTOS with the old correcting weight and SAm for the W — ^ fJV{y) 
decay distributions. 



of effects from complete matrix element to events generated with default PHOTOS. We can see 
also that distributions obtained from kernel of ref. [3 8 J and Y^x.e 1^/1^ close to each other but 
nonetheless distinct. Complete implementation requires control of the y* spin. The Bom level 
distribution for y* decay has zero at cos0^± = ±1. Close to these directions internal weight 
of PHOTOS necessary for exact matrix element becomes large, see figs. 1131141 Zero which is 
also present in distribution of Bom level W decay is of no such consequences for the weight 
distribution. 



6 Summary 

In this paper we have studied matrix elements for the 7* — )■ TT+TT^y and W — )■ Ivy processes. 
We have observed that the expressions can be separated into gauge invariant parts. In both 
cases, the part consisting of eikonal factor multiplying the Born level spin amplitude separates. 
This part contributes to the infrared singularity. In the case of 7* — )■ n^n^y the remaining part 
does not contribute to coUinear singularity but for — t- Zv it can be separated further: into 
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Figure 7: Ratios of results from PHOTOS with the excact correcting weight and SANC /or the W — t- ^v(y) 
decay distributions. 



the part proportional to the W charge which does not lead to any logarithmic contribution after 
integration, and the part proportional to the lepton charge, which do contribute to coUinear 
singularity and is identical to an analogous part for the Z — )■ /+/^ process (see eg. second 
or third part of formula 5 in [IZTl ). This is exactly the factorization property needed for the 
iterative solution used in PHOTOS to be valid for multiphoton emissions and processes discussed 
in present paper. 

For 7* — )• TT+TT^y, the factor, identified as Born level amplitude is not unique. As expected, 
ambiguity is proportional to photon momentum, and disappears in the soft limit. As a con- 
sequence of the ambiguity options were discussed. Their difference is of no practical conse- 
quences for the present work, but should be kept in mind if the obtained spin amplitude parts 
would be used as building bricks for amplitudes of more elaborated processes. 

We have identified dominant parts of spin amplitudes and used them for tests and for kernels 
of PHOTOS Monte Carlo. Exact matrix element was used for that purpose. We have found, that 
the whole matrix element for the process W — >■ /vycan be incorporated into the photon emission 
kernel. However for y* — Tt^Ti^y it is not technically straightforward, because of large weight 
events. The responsible term was identified and it may be interesting to point out that it is similar 
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Figure 8: Cosine of the angle of photon momentum with respect to spin axis in the decay W — )• fJv{y)- 
W boson is transversely polarized, = — 1- The black line corresponds to SANC results, while the red 
line corresponds to predictions o/PHOTOS with standard kernel. 
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Figure 9: Cosine of the angle of photon momentum with respect to spin axis in the decay W — t- /^v(y). 
W boson is longitudinally polarized, Jf' = 0. The black line corresponds to SANC results, while the red 
line corresponds to predictions o/PHOTOS with standard kernel. 



to the one obtained in a different calculation (Ref. [Tf| equation (68)). There, such terms were 
interpreted as contributing to running of QCD coupling constant. In present calculation the term 
involves complete kinematics of y* — t- tt+tt^y process. 
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Figure 10: Cosine of the angle ofmuon momentum with respect to spin axis in the decay W — )• /^v(y) is 
shown. W boson is transversally polarized, = — 1- The black line corresponds to SANC results, while 
the red line corresponds to predictions o/PHOTOS with standard kernel. Comparison of the distributions 
show that lack of spin effects in PHOTOS standard kernel results in up to 14% surplus of events for phase- 
space region close to cos 6^ = — 1. 
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Figure 1 1 : Cosine of the angle of the muon momentum with respect to spin axis in the decay W — )• /uv(y) 
is shown. W boson is longitudinally polarized, = 0. The black line corresponds to SANC results, while 
the red line corresponds to predictions o/PHOTOS with standard kernel. Comparison of the distributions 
show that lack of spin effects in PHOTOS standard kernel results in up to 4% surplus of events for phase- 
space region close to cos 6^ = ±1. 
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Figure 12: Cosine of the angle of photon momentum with respect to spin axis in the decay W — t- ^v(y) 
is shown. W boson is transversely polarized, = — 1, along this axis. The black line corresponds to 
SANC results, while the red line corresponds to predictions o/ PHOTOS with exact matrix element. 



Figure 13: Angular distributions for e^e^ — )• 71+71^(y) at 2 GeV center of mass energy. Results from 
PHOTOS with matrix element taken from l[38\l are given in red colour Matrix element Y,X.e l^/P used 
for results with green line. Logarithmic scale is used, but for the ratio (black line) linear scale is used 
instead. Fraction of presented events (i.e. with photons above 50 MeV) is respectively 4.2279 it 0.0021 
% and 4.2269 it 0.0021% of the total samples for the two programs. 
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(a) cosine of the photon angle with respect to beam, (b) cosine of the 7t+ angle with respect to the same 

charge beam. 
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Figure 14: Angular distributions for e^e^ — )• K^K^{y) at 2 GeV center of mass energy. Results from 
PHOTOS with matrix element taken from / liSl/ are given in red colour Y.x.e I^/'P used for results with 
green line. Logarithmic scale is used, but for the ratio (black line) linear scale is used instead. Fraction 
of presented events (i.e. with photons above 50 MeV) is respectively 4.2279 it 0.0021 % and 4.2271 it 
0.0021% of the total samples for the two programs. 
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(a) cosine of the photon angle with respect to beam is (b) cosine of the 7t+ angle with respect to the same 
shown. charge beam is shown. 



Any further investigation of the analogy, would require amplitudes of higher orders. Such 
effort can not be justified for scalar QED. From low energy point of view, terms like ([72)) 
or (fT4)) should be understood as genuine scalar QED process, and not part of real photon 
bremsstrahlung. Scalar QED is not supposed to be valid in the regions of phase space where 
these terms contribute significantly. It is of phenomenological interest to check this limit of 
scalar QED predictions by direct comparisons with data. Thanks to present work, higher order 
genuine bremsstrahlung effects for 7* — TT+TT^ycan be simulated with the help of PHOTOS and 
one can concentrate on confronting the data with these non-bremsstrahlung parts (fT2l) or (fT4l) . 

We have neither discussed here the interference with the photons originating from incoming 
beams, nor the interference between two consecutive emissions from the same charged line. 
The first effect, requires simultaneous treatment of initial-state and final-state bremsstrahlung. 
This is out of scope of work on PHOTOS alone, but spin amplitudes are already prepared. For 
discussion of interference of two emissions from the same charged line (and resulting uncer- 
tainties) second order matrix element is needed. Fortunately the structure of spin amplitudes 
for 7* — > TT+TT^y and W ^ Iv matches that of Z — )■ l^l^y [i33i . At present, we can only expect 
that these results on Z — )■ /+/ y in combination with our algorithm for matching consecutive 
emissions, hold for our y* -> n^n^ processes. Results of Ref. [ , 54j point that this expectation is 
well founded. 

We have not discussed virtual corrections. We assume, following Kinoshita-Lee-Nauenberg 
theorem ||481|49jl, that the dominant part can be included in a factor multiplying Bom amplitude 
and the correction to the total rate is free of any large logarithm. We leave this point for future 
work. 
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Finally, this paper provides numerical tests of PHOTOS Monte Carlo, in particular construc- 
tion for decays where Bom level cross section has a zero. This is of practical interest for users 
of the program and also a necessary step before any attempt of an extension to QCD. 



A Matrix element for W decay 

The matrix element of the process W^{Q,'k) — )■ /(/>/, ?1/)v(/>vAv)y(^,<5) has the form 



^' b(y{k,pi) - {ba{k,pi) +ba{Kpv)) 



2k- pi 



2k-Q 



Qi 



2k^, E Ul,{pi,mi,k,0,k,0)B^p,^X.{k,Q,p,) 
Qw 



2k. Q^^ (^t ,-P ' 2' ^)^"p,-;.v 0' 0' i^v, 0) 
+ t/jj (pi, m, , k, 0, k, 0)Bl^_^^ {k, 2, py) 



(32) 



where we use the following notation : 



2V2 



i/(piAi)ei^(e)(l+Y5)vb2,^2) 



Ux^,X2{Pumi,k,0,p2,m2) = u{pi,Xi)£y{k)u{p2,X2) 
hiX2^a{k,p) = U^^j^^{p,m,k,0,p,m), 



(33) 



Qi and Qw are respectively the electric charges of the fermion I and the W boson, in units of 
the positron charge, e!^ (^) and £^^^(2) denote respectively the polarization vectors of the photon 
and the W boson. An expression of the function U^^ in terms of the massless spinors and 
other notations can be found in [|52|. It is easy to check that the three components of the sum 
contributing to (|32l) are individually gauge invariant. Note, that the first component coincides 
with the amplitude in the eikonal approximation. 
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